A geometric model to diffraction of light by topological defects in liquid crystals 
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We study the light scattering by a hedgehog-like and linear disclination topological defects in a 
nematic liquid crystal. Propagating near such defects, the light feels an effective metric equivalent to 
the spatial part of the global monopole and cosmic string geometries. We obtain the diffraction pat- 
tern in terms of the characteristics of the liquid crystal, the scattering amplitude and the differential 
and total scattering cross section for the case of hedgehog defect. Studying the disclination case, a 
cylindrical partial wave method is developed. As an application of the previous developments, we 
also examine the temperature influence on the localization of the diffraction patterns. 
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I. INTRODUCTION 

Recently, a geometrical model for light propagation in 
nematic liquid crystals has been put forward by one of us 
and a collaborator In this model Fermat's principle 
is equated to the variational principle that determines 
geodesies in Riemannian geometry. This way, the path 
described by the light is interpreted as a geodesic in an 
effective geometry. Although one cannot gather more 
information than that which is already provided by Fer- 
mat's principle the tools of Riemannian geometry come 
in very handy to produce results that otherwise would 
be very complicated to get 0]. In this article we make 
use of the effective geometrical background in order to 
study diffraction of light by topological defects in nemat- 
ics. We look at a disclination and a point defect (hedge- 
hog), each of topological charge -1-1. The choice is due 
to their azimuthal symmetry which much simplifies the 
calculations. 

Topological defects in nematics, as in many physical 
systems, appear naturally in a symmetry-breaking phase 
transition. The high temperature isotropic phase has 
complete 3D rotational symmetry which is reduced to 2D 
rotational symmetry around the director axis of the ne- 
matic phase. Quenching the isotropic phase leads to the 
formation of domains with different orientations in the 
nematic phase and, consequently, to the creation of de- 
fects. Among these, the commonest are the disclinations 
and point defects named hedgehogs. Something analo- 
gous apparently happened at the beginning of the uni- 
verse: as the universe expanded and cooled down it went 
through a sequence of symmetry-breaking phase tran- 
sitions that left behind topological defects like cosmic 
strings and monopoles. These have been shown [4] to be 
analogous, from the point of view of light propagation, 
to disclinations and hedgehogs. As noticed by Bowick 
et al. 3, the physical principle that conducts the cre- 
ation of these defects in liquid crystals is the same found 
in the creation of cosmic strings for the universe: the 
mechanism of Kibble . 

There has been much effort put on analogous gravi- 
tational systems like the one described above in recent 



years [g, |7[ . Because these systems present similar math- 
ematical description, analogous models suggest the pos- 
sibility of reproducing selected gravitational systems in 
the condensed matter laboratory, allowing the study of 
aspects of general relativity not accessible to the cur- 
rent technology. For example, when a fluid flow has a 
downstream current, any phonon propagating in this sur- 
roundings feels a kind of a gravitational pull and black- 
hole-like equations rule the phonon's path Q]. In this 
example, the phonon propagation is described by an ef- 
fective Riemannian geometry, such that null geodesies of 
this geometry are traced by the phonons. Similarly, the 
use of an effective metric, as described in pj, will be done 
in the present article to study new, wave-like, aspects of 
the propagation of the light. 

In this article we study the propagation of light around 
topological defects in nematics from the point of view 
of wave interference and diffraction. We examine what 
happens to a plane wave which collides with a disclina- 
tion and on a hedgehog defect. We calculate the scatter- 
ing amplitude for both cases and obtain the respective 
diffraction spectra. 

This article is divided as follows. In the first part, 
we present the geometric model and the effective metrics 
used throughout this work. In the second and the third 
parts, we apply the partial wave method to the metrics 
generated, respectively, by hedgehog and disclination de- 
fects. On the hedgehog case, we obtain the diffraction 
pattern, its localization, the scattering amplitude and 
the differential and total scattering cross section, while 
on the disclination case, we emphasize on the localiza- 
tion of the diffraction pattern. In the fourth part, we 
analyze the temperature influence on the localization of 
the diffraction patterns and, in the fifth one, we present 
our conclusions. 



II. GEOMETRIC MODEL 

Observing the light propagating between two isotropic 
media, no effect in its path is noticed, besides those fore- 
seen by Fresnel's equations. However, considering the 
passage of a light ray from an isotropic medium to an 
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anisotropic one, for example, a liquid crystal in the uni- 
axial nematic phase (where the uniaxial attribute is rep- 
resented by the director vector n), one observes a split 
from the original ray to two other rays with perpendic- 
ular polarizations among themselves: the ordinary and 
extraordinary rays. The former's electric field polariza- 
tion is perpendicular to n, as long as the latter's one has 
a component lying on n. To calculate the propagation 
of a monochromatic wave, through Maxwell's equation 
in a medium without currents and charges, one notices 
i that the refractive index for the extraordinary ray is 
given by 

N^{(3) = nl cos^ (3 + nl sin^ f3, (1) 

where /3 is the angle between the director n and the 
Poynting vector S of the extraordinary ray (that is the 
actual direction of light propagation because S indicates 
the energy flow), and rig and Uo are, respectively, the 
refractive indexes of the rod molecules which constitute 
the liquid crystal. 

With this refractive index, it is possible to calculate 
the path of the extraordinary light applying the Fer- 
mat's principle. This says that the path traveled by light, 
among two points A and i?, will be that which minimizes 
the integral 

f Nedl, (2) 

J A 

that possesses dimension of length, using a parameter I 
along the path. Examining Riemannian geometry, the 
light's trajectory among two points will be that which 
minimizes the line element ds given by 

3 

ds^ — gijdx^dx\ (3) 

where, — {r, 0, in spherical coordinates, and gij 

are the components of metric tensor g. Thus, it is possi- 
ble to derive the geometric model by noticing that both 
Fermat's principle and the calculus of geodesies in Rie- 
mannian manifolds are questions of minimizing physical 
amounts. Bearing this in mind, if we compare (Ngdl)'^, 
from Fermat's principle, with ds^, from a generic Rie- 
mannian manifold, we interpret the path traveled by light 
as a geodesic in an effective metric, that is, we are led to 
the following identity 

Nldl^ = ^g,^jdx'dx^ . (4) 

This way, it is possible to calculate the effective metric 
felt by the light finding cos (3 and sin /3, in ([1]) and sub- 
stituting ll]) in ([4]). This procedure is done in [l| to the 
hedgehog-like defect and disclinations in the Inematic 
phase of a liquid crystal. The result is that an effective 
metric felt by light in the presence of a hedgehog-like 
defect is 

ds^ = dr^ + {dO^ + sin^ 6 dtp"^) , (5) 



where b = and in the presence ofa{fc = l,c = 0} 
disclination, 

ds^ = dp^ + b^p^d(j? + dz"^. (6) 

The constants k and c relate the angular polar coordi- 
nate (f) with the angle (f between the major axis of the 
liquid crystal molecule and the x-axis Equation ([5|) 
is identical to the spatial part of the global monopole's 
line element, being h the deficit of solid angle. Equation 
^ , as well, is identical to the spatial part of the cosmic 
string's line element, where b is the deficit of polar angle 
of the cosmic string. With this result in hands, the cal- 
culation of the light diffraction in the presence of both 
liquid crystal defects will be developed in the next two 
sections. 



III. HEDGEHOG DEFECT 

In this section, we will apply the method of partial 
waves ^0|] to calculate the light diffraction inside a liquid 
crystal in the nematic phase, being this altered by the 
presence of a hedgehog-like defect. For simplicity, we 
will consider the light as a scalar wave, allowing the use 
the d'Alembert's wave equation: 

v.v^$ = (y^ff'^a,<f) = 0, (7) 

where gij are the elements of the metric represented by 
([5|) and g is the determinant of gij. However, a comment 
must be done. The geometric model used here generates 
only the spatial part of a pseudo-Riemannian metric, i.e., 
for us, time is a parameter. This fact can mathematically 
be represented by the addition in © of the element —dt^. 
With this correction, the metric that we will apply in ([7]) 
will be 

ds^ ^ -~dt^ + dr^ + b^r^ (de^ + sin^ 9 dcj)^) . (8) 
Considering light with a given frequency uj, 

<i>(t,r,0,0)=e-^'-*^(r,0,^), (9) 
and using dH) and dH in ([7]), we obtain 

(A - c^2) 7/^ = 0, (10) 

being A written by 

A = — H , 

and the operator L is recognized as the angular part of 
the Laplacian applied to Lorentzian space (flat space) in 
spherical coordinates. At this point, we can use separa- 
tion of variables at spherical harmonics over ■0(r, 0, 0). 
However, the spherical symmetry implies no dependence 
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in the (j) coordinate (m = for spherical harmonics). 
Thus, the spatial part for light wave equation is 



(11) 



1=0 



As the main quantity of the method of partial waves is the 
phase shift, the calculation of the term ai in the previous 
equation will be ignored. Thus, when applying (jlip in 
(fTO|) . we will have the following radial equation: 



1(1 + 1) 



And after the change of variables Ri — -^j^ in the previ- 
ous equation, we find: 

rG 
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The solutions of this equation can be described by first 
kind Bcsscl functions [ll| 

Gi = Jn{i){(^r), 

where the phase of this function is 

1 

2 
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We observe that when b = ne/uo — 1 (see Eq. ([S])) we 
get n{l) — I + ^, representing a wave propagating in flat 
space, i.e., in the absence of the defect. Comparing the 
phases of the wave functions in the presence of ^ ^"■^'■^ ^ 

and in the absence of the defect (^ (^l + i)) , we find the 
following phase shift 
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(12) 



From the partial waves method [10|, we use the expres- 
sion to scattering amplitude, 

^ oo 

/(^) = ?- E(2^ + - micosO). (13) 

I 



For a given numerical value to b at (jl2p and substitut- 
ing into (fT3|) . the latter equation allows us to calculate 
numerically the differential scattering cross section a (6) 
(FiglT]), i. e., the angular distribution of scattered light, 
by the following expression 

aie,<p)^\fm'- (14) 
The azimuthal symmetry of the scattering generates a 
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FIG. 1: Scattering amplitude for 6 = 0, 9, = 1 and equation 
(fT3l) truncated at Z = 600. 



annular diffraction pattern, with its maximum located at 
the angle indicated by Fig. [1] This diffraction ring is the 
main characteristic of a scattering in the proximity of a 
global monopole. 

Using the method presented in [l^], on the b^ ^ 1 
regime , we can develop an analytical expression for the 
angle that locates the diffraction ring and for the total 
scattering cross section crtot given by the optical theorem. 
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Initially, making ( — I + ^ and 



(15) 
in (HID, we 



have 
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We can expand this last equation in terms of when 
considering 6^ ~ resulting 



(16) 



As is small, the two first terms of the previous expan- 
sion are enough relevant for approximate calculations. 
We can apply (dH) in ^ and find the zeroth, f'^°^6), 
and first, f'^^\0), order terms for the factor a^, for the 
scattering amplitude. These two terms can be written 
through the function 



h{9,a) ^^e'^*"^^') Pi (cos 6*) = 



(=0 



•\/2(cos7rQ! — cos6') 
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related with the generating function for the Legendre 
polynomials where a = 1 — -j. The first term of the 
expansion of the scattering amplitude is proportional to 
the derivative of h{9, a) in relation to a: 
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(17) 



Substituting this result in the equation (llSp . we find the 



first term of the expansion atot — c 



(0) 



r(l) 
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fc2 sin^ ^ 



(18) 



The second term of the scattering amplitude, after con- 
sidering the order O(a^) in (fT6)) . is proportional h{9,a): 



/(i)(0) 



1 



26A; Y^2(cos ira - cos( 



(19) 



And, again, applying this result in the optical theorem, 
we find: 



.(1) 



1 



sin ^ ■ 



(20) 



As it can be seen, equations ([T7| and (fT9|) have singular 
behaviors for the angle 9q = ira. Applying this result in 
p4p . we find an annular diffraction pattern located in 
00 = 7r(l — 1), valid only for 6^ w 1. However, when 
comparing this localization with the one gotten by means 
of numerical calculation developed by a direct application 
of p2p in (|13p . we see that both predict the same position 
to the diffraction ring for a given b. This means that the 
terms f''^\0) and /'^■'(0) are indeed predominant in the 
calculation of the accurate value of f{9). 



IV. DISCLINATION 

As we deal with a cylindrically symmetric problem, 
i.e., V{r) = V{r), then becomes natural, for the study 
of the light diffraction, the application of the method of 
cylindrical partial waves, which will be exposed as follow. 

The method of the cylindrical partial waves will be de- 
veloped similar to [l0| for the spherical case. The wave 
problem in this section is related to a defect that is pro- 
longed infinitely in the direction z, allowing only two 
space variables (p, <j)) . Consequently, the wave equation 
will be given by 



(A 
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0, 



(21) 



X, e*'^^, and a scattered wave (being ignored the prob- 
lem of the normalization). In this way, the behavior of 



V 



(diff) 



(r) for great values of p will be 

,f//)(r)^e'^- + A.(0)^ 



(22) 



In order to apply the method of partial waves, we need 
the solution of equation (|2ip in the absence potential, 
that is, the free cylindrical wave ip^^j^p, 0), 



lim ^f]{p,(f)) 
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l=cos(fcp-|(^-i,,x 



X e 
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(23) 



the solution of the equation (PT|) in the presence of a 
symmetrical potential V{v) = V{p), i.e., the cylindrical 
partial waves ipk.iip, (j^), 



lim %l;k,i{p,(t)) 

p— i-OO 



being 5i the phase shift regarding the free wave (|23p . and 
the plane wave expanded in free cylindrical waves using 
the Jacobi- Anger's expansion 



^ikx ^ikp cos 



lij, 



(25) 
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where eo = 1 and e; = (l + e^^"-*) (I > i). 

With these three elements, it is possible to expand 

v'f'^^^^ (r) in a series of cylindrical partial waves in such 
a way that the coefficients of (P5|) will be used so that 



■"k 



(r) reproduces a plane wave in p — > oo: 



„idiff) 
"k 



(r) = ^cii^kAP^' 



1=0 



(26) 



For p oo and using the exponential form of the func- 
tion cosine, the substitution ([M)) in results in 



where the factor p~ guarantees that the total flow of en- 
ergy passing through a circle of radius p is independent 
of p for great values. Like this, the wave function that 
represents the scattering state, v^i^^^^^^^\ will be the com- 
position of the plane wave, for instance, spreading along 
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Multiplying the right side of the last equation by the 
factor without physical meaning e*^' and using e^**' = 
1 + 216"^^' sin (5;, one attains 

oo 

/=o 

^,(fcp-^(;-i))^.^i^, Sin(5; \ 

^ j ■ 

Recognizing the first term between parentheses as be- 
ing the asymptotic form of (|25|) and evidencing the term 
^ly^ in the second term, the comparison of this with (1^^ 
identifies fki4>) as being 

oo 

fk{4>) = -=^ieT^qe^'*sin5ze^'5'. (27) 

Applying (|77|) in (fH)) we have the differential scattering 
cross section for cylindrical defects with angular symme- 
try. 

Now we apply the theory developed so far to calcu- 
late the scattering amplitude and the differential scatter- 
ing cross section for the scattered light by a disclination 
with effective metric represented by ([5]). This proceeding 
seems to be a simple case of the hedgehog defect seen 
from its equatorial section (0 = 0), that can be verified 
comparing the equations (O and In this manner, our 
starting point will be again the application of the effec- 
tive metric for the disclination ([6]), with the term —dt^ 
correcting the time part, in the wave equation ([7]). 

Due to the symmetry in the direction z, it will just 
study the propagation along the coordinates p and (a 
different approach of the same problem of the cosmic 
string can be seen in [l3|)- Therefore, using a separation 
of variables similar to we found the following wave 
equation in the presence of a disclination 

where it will be used ipiPT4') = YldLfi^i^iip)^^^'^ ■ With 
the latter expansion, the equation to the radial coordi- 
nate will be 

\ or^ r or r'^ ) 

being i?.(p) = Ji,^{p),iyi — ^. Using this index, one can 
find the phase shift on the way to equation (f24|) and, 
consequently, the scattering amplitude (I27p and the dif- 
ferential scattering cross section (ITU) . 

The phase shift of the scattered wave by the disclina- 
tion © is found to he Si = ^ (l - i)- Fig. [U shows 
the differential scattering cross section with the same pa- 
rameters of the Fig [1] Analogous to the latter. Fig. [2] 
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FIG. 2: Scattering amplitude for 6 = 0, 9, A; = 1 and equation 
((27|) truncated at Z = 600. 

indicates, for the disclination case, a singular behavior of 
around 0o = Tra, where a = (I — ■^), generating a 
luminous line representing the diffraction pattern about 
Like this, one could imagine the disclination as a line 
of punctual hedgehog-like defects, or from the gravita- 
tional point of view, a cosmic string as a line of global 
monopoles. 



V. TEMPERATURE DEPENDENCE 

It is well known the temperature dependence that the 
refraction indexes of a liquid crystal molecule possess fT^, 
[l6j . We will use this fact to calculate the influence of the 
temperature in the position of the diffraction pattern for 
the previously studied defects. 

As exposed in [T^, |T3], the refractive indexes for or- 
dinary and extraordinary rays traveling in a liquid crys- 
tal molecule can be treated by a four-parameter model. 
Thus, the refraction indexes are expressed by 

ne(r) « A-i?T+^^(^l--j (28) 
n„(T) « A-BT-^-^(\-^ , (29) 

where A is a dimensionless constant, i? is a constant with 
(K^^^ unit, (An)^ is the birefringence of the material 
in the crystalline state, /? is a material constant and Tc 
is the transition temperature between the isotropic and 
nematic phase. Using this model for 5CB (4-cyano-4- 
n-pentylbiphcnyl) , that has = 306.6K, and a light 
beam with wavelength A = 589 nm, one obtains A = 
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FIG. 3: Temperature dependence for the position of diffrac- 
tion pattern, 0o = tt (l — i) , for 5CB. The dots are numerical 
results. 

1.7674,5 = 5.79 x 10-^K-\ (Aji)^ = 0.3505 and /3 = 
0.1889 O]. Applying ([28]) and dSH]) in jS]) and dSl), one 
obtains a new expression for b and, consequently, the 
temperature dependence of the diffraction pattern of a 
hedgehog-like and discUnation defect by the relation (j>Q = 

Fig. [3] shows the behavior of the diffraction pattern 
position for 5CB in its nematic phase (295 - 306.6 K 0) 
regarding the two previously studied defects. One ob- 



serves the low sensibility on the position at the beginning 
of temperature range (« 295 K) and a high sensibility at 
the end of the temperature range (« 306 K). 



VI. CONCLUSION 

The relationship between the Fermat's principle and 
its description by an effective metric allows us to acquire 
information about the optical properties of the liquid- 
crystalline material calculating the position of the diffrac- 
tion pattern resultant from the light scattering over a 
hedgehog-like or disclination defect in the nematic phase, 
because 9o = 7r(l — i), where b = With these 

two relations, it is possible to measure me temperature 
by the diffraction pattern localization of the light scat- 
tered by hedgehog-like or disclination defects, mainly 
near the high temperature end of nematic phase. Thus, 
the present work serves as example of applicability of 
analogue models, in special, using liquid-crystalline ma- 
terials. 
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